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Abstract. An abstract should be given 



The He I A5876A absorption line recently discovered in 
the spectrum of the companion to the millisecond pulsar 
PSR J1740-5340 is tentatively attributed to electron im- 
pact excitations due to the irradiation of its atmosphere by 
7-rays emitted by the pulsar's magnetosphere. Numerical 
calculations, similar to those carried out previously for 
Type lb SNe, indicate that a pulsar beam with photon 
energies ~ 1 MeV gives rise to a A5876A line of the ob- 
served strength if the beam's spin-down conversion effi- 
ciency approaches 1%. However, a significant difficulty for 
the proposed mechanism is the strength of the singlet line 
at A6678A. Compared to the corresponding triplets, sin- 
glet lines are weak because of the loss of excitation when 
photons emitted in the series n^P — > l 1 ^ ionize hydrogen 
atoms, an effect absent in the hydrogen-free atmospheres 
of Type lb SNe. 

Key words, pulsars: individual (PSR J1740-5340) - radia- 
tive transfer - stars: atmospheres - line: formation 



1. Introduction 

The binary millisecond pulsar PSR J1740-5340 was re- 
cently discovered in the globular cluster NGC 6397 by 
D'Amico et al. (2001a) with the Parkes radio telescope. 
They reported a spin period of 3.65ms, a 1.35 day circu- 
lar orbit, and eclipses at 1.4 GHz that last for more than 
40% of the orbital period. Subsequent Parkes observations 
(D'Amico et al. 2001b) allowed an accurate timing solu- 
tion yielding the spin-down rate as well as the binary's 
coordinates. 

With the position known, Ferraro et al. (2001) used 
archived images to obtain a secure optical identification 
with a star (WF4-1) previously listed as a BY Dra candi- 
date (Taylor et al. 2001). The pulsar's companion (named 
COM J1740-5340 by Ferraro et al. 2003) has a luminos- 
ity close to that of the cluster's main-sequence turnoff 
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stars but its colour is anomalously red. In addition, the 
companion is variable with a range of ~ 0.2mag. This 
variability has the same period as the binary orbit and 
matches that due to the tidal distortion of a star filling its 
Roche lobe. Ferraro et al. (2001) emphasize that, in this 
regard, PSR J1740-5340 differs from two other eclipsing 
millisecond pulsars, where strong photometric variability 
is clearly due to the heating of the companions' facing 
hemispheres by the pulsars' emission. The absence of this 
expected heating effect is also reported by Kaluzny et al. 
(2003) and by Orosz & van Kerkwijk (2003). 

Although not apparent in broad-band light curves, the 
effect of the pulsar's emission on its companion has almost 
certainly been detected spectroscopically by Ferraro et al. 
(2003). From observations made to determine the spec- 
troscopic orbit, these authors made the serendipitous dis- 
covery of He I absorption lines at 5876 and 6678A. Given 
that the companion's effective temperature T e // — 5530K 
(Ferraro et al. 2003) is substantially less than the T e ff 
(<; 10, 000K) at which the He I lines are present for main 
sequence stars, this discovery is a major surprise and ob- 
viously significant. 

Ferraro et al. (2003) explain the presence of He I 
absorption in terms of a narrow strip across the illu- 
minated face of the companion within which the local 
T e ff ~ 10,000K due to heating by beamed pulsar radia- 
tion. But outside this strip, irradiation is negligible and so 
the local T e / / is then ~ 5000if . By a suitable choice of pa- 
rameters, Ferraro et al. imply that the observed strengths 
of the A5876A line can reconciled with the absence of a 
photometric heating effect. In a very recent paper, this 
group (Sabbi et al. 2003) have published more of their 
spectroscopy and repeated their conclusion that the pul- 
sar's emission pattern must be highly anisotropic. 

Although the model proposed by Ferraro et al. is not 
implausible, another effect suggests itself if we note that 
COM J1740-5340 has two points in common with Type 
lb SNe. The first is that the spectra of these SNe are also 
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remarkable in exhibiting He I absorption lines despite cool 
photosphcric continua (Harkness et al. 1987). 

The second point in common - the presence of high- 
energy photons - follows from the explanation for the he- 
lium lines in the Ib's. In the standard interpretation (Lucy 
1991, hereafter Paper I; Swartz et al. 1993 ), a huge over- 
population of the n — 2 levels is caused by impact exci- 
tations and ionizations of helium atoms by non-thermal 
electrons. These energetic electrons are the cascade prod- 
ucts created by the multiple Compton scatterings of 7-rays 
emitted by the radioactive nuclei 56 Co. There is a signifi- 
cant energy density of such 7-rays in the reversing layers 
of these SNe because of mass loss preceding the explosion 
and the mixing out of 56 Ni following the explosion. 

In view of these similarities, an obvious conjecture is 
that 7-rays emitted by the magnetosphere of PSR J1740- 
5340, or arising when the pulsar's relativistic wind en- 
counters the companion's wind, irradiate the atmosphere 
of its companion, overpopulate the n = 2 levels of He , 
and hence give rise to detectable absorption lines. The in- 
vestigation of this conjecture is the subject of this paper. 

2. An irradiated model atmosphere 

In this section, the problem is defined and then the steps 
needed to calculate the He line spectrum are briefly de- 
scribed. The discussion follows closely that given in Paper 
I for Type lb SNe. 

2.1. The problem 

In investigating the non-thermal excitation of helium in 
the atmosphere of COM J1740-5340, we face several ma- 
jor uncertainties. First and foremost, there is no direct 
evidence that PSR J1740-5340 emits a powerful beam of 
hard radiation. Although the object is detected in X-rays 
(Grindlay et al. 2001, 2002), its X-ray luminosity Lx is 
weaker by a factor ^17 than expected for its spin-down 
luminosity E from the correlation Lx ~ 10~ 3 -E found for 
millisecond pulsars in the field (Becker & Triimper 1997). 
However, the detected X-rays are probably not emitted 
from the pulsar's magnetosphere but rather from the in- 
teraction region between the pulsar's relativistic wind and 
the mass loss from the companion (Grindlay et al. 2002). 

Despite this absence of evidence for a beam of hard 
radiation, we assume that such a beam exists and that 
it sweeps over at least part of the companion's disk. 
Accordingly, the problem investigated in this paper is to 
study the formation of the Hei spectrum at a typical point 
on the companion's disk lying within the strip swept out 
by the beam. 

The problem as posed refers to irradiation by a pen- 
cil beam emitted by the pulsar's magnetosphere. But the 
calculation is also relevant to irradiation by an extended 
source of hard radiation, such as that due to the colliding 
winds of the pulsar and its companion. The single beam 
then represents a one-point quadrature approximation for 
the effect of the extended source. 



2.2. The incident beam 

Because everything about the beam is conjectural, we take 
it to be mono-energetic, comprising photons of single en- 
ergy eoj the beam's effective photon energy. If the beam's 
flux as a function e were known, eq could be defined as the 
value of e such that the mono-energetic beam deposits the 
same fraction of its flux in the reversing layer as would the 
actual beam. Nevertheless, since the vertical distribution 
of the deposited energy would not be exact, the concept 
of an effective photon energy is not rigorous. 

The phase-averaged flux incident on the companion's 
atmosphere is conveniently written as 



To = no x-q x 



E 



47T 



(i) 



where [i is the incident direction cosine, a is the orbital 
separation, 77 is the conversion efficiency of E into beamed 
radiation, and fi is the solid angle swept out by the beam. 

The adopted numerical values are as follows: E — 
1.4 x 10 35 erg s^ 1 (D'Amico et al. 2001b); a = 4.33 x 10 11 
cm , corresponding to component masses M\ — 1.5 Mq, 
M 2 = 0.26 Mq (Ferraro et al. 2003); and ft = 1 stera- 
dian. This last choice is standard practice in view of our 
ignorance concerning beaming fractions. 

Inserting these values, we have T a = 7.48 x lO 11 /^?? 
erg cm~ 2 s _1 . 

2.3. Model atmosphere 

A full NLTE investigation of the response of the com- 
panion's atmosphere to irradiation by a beam of 7 rays 
is not justified in view of the above-mentioned uncer- 
tainties. Accordingly, the approach adopted is to em- 
bed a NLTE calculation for He atoms in an atmo- 
sphere that is otherwise unperturbed by irradiation. 
Moreover, a severely simplified model is adopted for this 
unperturbed atmosphere. Specifically, a plane-parallel, 
Schuster-Schwarzschild model is assumed comprising a 
black-body emitting photosphere at Rosseland mean opti- 
cal depth r — 2/3 and an isothermal reversing layer with 
temperature T = 0.93T e //. On the assumption of hydro- 
static equilibrium and constant Rosseland mean opacity 
k,r, the variation of gas pressure in the reversing layer is 
given by P = t x <?/k_r, where g is the surface gravity. 
The adopted constant opacity hr is the solution of the 
equation P = 2/3 x g/n R {P,T e ff). 

The parameters of the model atmosphere are T e ff — 
5530K (Ferraro et al. 2003), logg = 3.44, X=0.75 and 
Y=0.25. The value of logg corresponds to radius R 2 = 
l.6R Q (Orosz & van Kerkwijk 2003) and mass M 2 = 0.26 
A4q. The adopted chemical composition neglects metals 
since [Fe/H]~ —2 (Ferraro et al.) and assumes that enve- 
lope stripping has not yet exposed He-enriched zones. 

With these parameters, the equation for the opacity 
has the solution = —1.13 dcx. This was obtained by 
iterative interpolation in Table 6 of Alexander & Ferguson 
(1994). Although this Table is computed for Z = 0.02, the 
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H~ ion is the dominant source of opacity for the given pho- 
tospheric parameters, with metals as minor contributors 
to the electron density n e . 

In the subsequent calculations, the reversing layer is 
discretized into 50 layers of width Alogr = 0.1, within 
each of which physical variables are constant. The result- 
ing small optical depths of the outermost layers ensure 
that line centres are accurately calculated. 

2.4. Beam loss rates 

The first step in calculating the He spectrum is to derive 
the rate at which the beam deposits energy in the reversing 
layer in the form of energetic electrons. Each incident 7- 
ray creates a cascade of non-thermal electrons by multiple 
Compton scatterings off free and bound electrons. 

At optical depth r, the rate at which 7-rays are losing 
energy by Compton scattering is 

C{t) = A-k J f e a e J e de (2) 

where a e is the Compton scattering coefficient, f e is the 
expected fraction of e transferred to the Compton electron, 
and J e (r) is the mean intensity at photon energy e. Note 
that here, and throughout the paper, rates refer to unit 
volume. 

As in Paper I, a Monte Carlo code is used to simulate 
this process. If the incident flux To at r = is represented 
by M identical photon packets containing 7-rays of energy 
eo and entering the atmosphere in time interval At, then 
each packet's initial energy wo is given by M x wo/At = 
To- These packets are followed as they propagate within 
the reversing layer until they escape at r — or penetrate 
below the photosphere at r = 2/3, where they are assumed 
to thermalize. 

When the scattering histories of all N packets are com- 
plete, the loss rate in the m-th layer is derived with the 
Monte Carlo estimator (Lucy 1999, 2003) 

Here V is the volume of a vertical column of unit cross 
section through the m-th layer, and As is the pathlength 
of a packet of energy w £ between consecutive events, where 
an event is either a Compton scattering or a crossing of a 
layer's boundary. The summation is over all pathlcngths 
in layer m. 

Monte Carlo convergence experiments show that this 
estimator gives fractional errors < 10~ 3 for Af = 10 6 , the 
sample size adopted here. Note that this estimator returns 
a non-zero value for the loss rate even if no 7-ray packet 
undergoes a Compton scattering in layer m. 

Because C is proportional to 77, it is convenient to com- 
pute its value C 1 for 77 = 1. The energy deposition rate in 
layer m is then 77 C] n (eo 1 Mo)- 



2.5. Deposition fractions 

The next step is to determine the rates of impact exci- 
tations and ionizations of He atoms. To achieve this, we 
must follow the slowing down of fast electrons as they in- 
teract with matter in the reversing layer and record the 
energy transfers into the various deposition channels. 

This calculation is greatly facilitated by two approx- 
imitions adopted in Paper I. First, noting that the stop- 
ping distance for an MeV electron is ~ 10 -2 that for an 
McV photon (Colgate, Pctschek & Kriese 1980), we as- 
sume that the electrons deposit all their energy in situ. 
The second approximation (Meyerott 1978) makes use of 
the asymptotic near constancy of Dk(E), the fraction of 
an electron's kinetic energy E that is deposited in channel 
k. Thus we need only compute the deposition fractions D k 
for £t = iokcV since D k {E) ~ d{ for E > E^ . 

The values D\ are obtained by solving the Spencer- 
Fano equation - Eq (2) in Paper I. This requires that we 
identify the stopping mechanisms and provide their cross 
sections. In Paper I, this equation was solved for a weakly- 
ionized pure helium gas, so that the stopping mechanisms 
were impact excitations and ionizations of He atoms and 
Coulomb scattering with free (thermal) electrons. Now the 
equation must be solved for a weakly-ionized H-He gas, so 
the stopping mechanisms to be added are impact excita- 
tions and ionizations of H atoms. 

The treatment of He impacts and Coulomb scattering 
remain as described in Sect. 2.3 of Paper I. The added 
impact cross sections for H are the approximate formulae 
of Johnson (1972), with the H atom represented by 14 
bound levels with principal quantum numbers n = 1 — 14. 
In addition, the probability distribution for the energy E s 
of the secondary electron when H is ionized is taken to be 
oc {E 2 S + J 2 )" 1 with J = 8eV (Shull 1979). 

With composition specified, the deposition fractions 
D\ are functions of just one variable x e , the fractional ion- 
ization of H. In conformity with the assumption (Sect. 2.3) 
that the atmosphere is unperturbed by irradiation, elec- 
tron densities are computed from Saha's equation. Then, 
neglecting the variation with height, we set x e = x* e = 
5.2 x 10~ 5 , its value at t = 1/3. With this approximation, 
the deposition of non-thermal energy throughout the re- 
versing layer is governed by the single vector D\.{x* e ). 

Summing the appropriate elements of this vector, wc 
find that 11.7% of the deposited energy is converted di- 
rectly into heat via Coulomb scattering, 82.0% is con- 
verted into ionization and excitation energy of H, and the 
remaining 6.3% is converted into ionization and excitation 
energy of He. 

2.6. Statistical equilibrium of helium 

The rate at which non-thermal energy is absorbed in excit- 
ing He atoms to level i in layer m is ijC]^ x D\ . Accordingly, 
if Wi is the energy of this level relative to the ground state, 
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the rate of impact excitation (ionization) of level i in layer 
m is 



V x £ m (e ,Mo) x 



(4) 



If AjiUj denotes the rate at which transitions j — » ?' 
occur due to conventional radiative and collisional (ther- 
mal) processes, then, with impacts included, the equations 
of statistical equilibrium for He levels i = 2, 3, . . . , k are 



(AjiUj - AijUi) 



(5) 



When the constraint Yli n i — n He is added, we have k 
equations in the k unknowns rij. This system must be 
solved in each layer of the atmosphere in order to com- 
pute the emergent profiles of the He lines. But since the 
radiative coefficients in Eq. (5) depend on the solutions of 
the line transfer problems, iterations are required (Sect. 
2.8). 

The data sources for the A coefficients in Eq. (5) are 
as stated in Paper I, except that the collision strengths of 
transitions between the lowest five levels of He have been 
updated using the calculations of Sawey & Berrington 
(1993). 

2. 7. Escape probabilities 

For most of the He lines, we can assume that line photons 
in the reversing layer interact only with their own transi- 
tion - i.e., continuum process can be neglected. But this 
is not true for transitions to the ground state since the 
emitted photons can ionize hydrogen and the atmosphere 
is optically thick in the Lyman continuum. 

For the adopted atomic model, this is relevant for the 
permitted ground state transitions n 1 P° — ► l 1 ^ at 584, 
537 and 522A. Since the optical depths for these transi- 
tions are also huge, we can assume that when such a line 
photon is emitted it is absorbed in situ cither by the line 
itself or by a H atom. Accordingly, we adopt an escape 
probability formalism for these three transitions and only 
consider thermal motions in computing Doppler broaden- 
ing. 

If <j> v denotes the normalized absorption profile, a line 
photon emitted with frequency v sees a line absorption 
coefficient l v = a v rij, where 



_ 7re^ 

®>v Jii 

m e c 



1 



9ji 



9i n ] 



(6) 



and a continuum absorption coefficient k v — aniy) nn, 
where we assume all H atoms are neutral and in the ground 
state and neglect the correction for stimulated emission. 
The probability that this particular photon will be lost to 
the line is therefore k v /(k v +t v ). Thus, if we average over 
the emission profile, the escape probability is 



Pa 



b v dv 



(7) 



Note that the emission profile is taken to be identical to 
the absorption profile on the assumption of complete re- 
distribution (Mihalas 1978, p.29). 

This treatment is applied to the permitted transitions 
i — > j for j = 1. As a result, the radiative contribution of 
these transitions to Eq. (5) simplifies to PijUi, where 
Aij is the Einstein coefficient for spontaneous emission. 

For the isothermally stratified reversing layer, these 
escape probabilities are independent of r. With <j> v given 
by the Voigt function assuming pure radiation damping, 
we obtain /3 y = 1.8 x 1CT 3 , 4.4 x 1CT 3 and 9.2 x 1CT 3 for 
the lines at 584, 537 and 522A, respectively. 

An escape probability treatment is also adopted for re- 
combinations to the ground state. The emitted photon will 
travel a negligible distance before it photoionizes a H or 
a He atom. Since the recombination photons are strongly 
concentrated to the He ionization threshold at v>t, the es- 
cape probability is accurately given by (Osterbrock 1974, 
p.23) 



&1 



n H a H {vT) 



n H a H (y T ) + riH e aHe(v T ) 



(8) 



Accordingly, if ol\ denotes the total recombination coeffi- 
cient for the ground state, it is replaced in the equations 
of statistical equilibrium by S3 K \a.\. With the He threshold 
cross section 7.48 x lCT 18 cm 2 of Fcrnley et al. (1987), we 
obtain /3 Kl = 0.68 

2.8. Line formation 

Subordinate transitions are not treated with escape prob- 
abilities but by solving the transfer equation. For such 
transitions, the radiative contribution to Eq. (5) is A^rii — 
(BjiUj — Bijm)Jij, where 



Jij = / J v 4> v dv 



(9) 



is the profile-averaged mean intensity. This quantity is 
obtained by numerical integration after solving the line 
transfer equation 



dz 



-l v {h - S) 



with source function 



2hv z 



- 1 



9ji 



(10) 



(11) 



for a grid of frequencies spanning the transition. Since high 
precision is not required, Eq. (10) is solved only at /i = 
±1/2 giving J„ = (7+ + I~)/2. The boundary conditions 
on the inward and outward beams are I~ = at r = 
and 1+ = B v {T e f f ) at r = 2/3. 

As in Sect. 2.7, the line profile </>„ is given by the Voigt 
function assuming pure radiation damping. But in view 
of the moderate optical depths of the subordinate transi- 
tions, a microturbulent contribution of 3 km s _1 is now 
included in the Doppler parameter. 
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A simple iterative procedure is adopted to obtain mu- 
tually consistent solutions of the equations of statistical 
equilibrium and line transfer. Eq (5) is first solved with 
Jjj = B v (T e f f)/2, its optically-thin limit. The resulting 
level populations are then used to obtain improved values 
of Jij from the transfer equation. This process is repeated 
until the average of the absolute fractional changes of the 
level populations in the atmosphere are < 1CP 6 . 

When the iterations have converged, each line's emer- 
gent flux is F„(0) = I„(0), from which its residual in- 
tensity r v — F v (0) I B u (T e f f) and equivalent width W\ = 
J(l — r v )d\ are readily computed. 

2.9. Hydrogen lines 

With only minor modifications, the code described above 
can also be used to compute the H line spectrum of the 
irradiated atmosphere of COM J1740-5340. The vector 
D\.(x*) already contains elements giving the fractions of 
rj £} m used to excite and ionize H atoms. Accordingly, 
with the standard assumption of in situ re-absorption of 
Lyman- continuum and line photons, the H line spectrum 
is similarly obtained by iteratively solving the statistical 
equilibrium and transfer equations of Sects. 2.6 and 2.8. 
One point of difference is that impact ionizations of H are 
allowed to increase n e from its LTE value. 

3. Numerical results 

In this section, the simple treatment of line formation in 
an irradiated atmosphere developed in Sect. 2 is applied 
to compute the strengths of He lines in the spectrum 
of COM J1740-5340. With atmospheric parameters fixed 
(Sect. 2.3), line strengths depend on /io, r\ and eo- Selecting 
a typical point on the disk, we set = 2/3 and thus limit 
the investigation to exploring the atmosphere's response 
to a beam with conversion efficiency rj and effective photon 
energy e . 

3.1. Deposition efficiency 

Line strengths depend on e only through the function 
£m( e o> Ho) defined in Sect. 2.4. From the Monte Carlo 
calculations made in determining this function, it is infor- 
mative to compute the fraction of incident flux deposited 
in the reversing layer, the fraction deposited below the 
photosphere at r = 2/3, and the fraction emerging back 
into space at r = 0. 

Fig. 1 plots these deposited fractions as a function of 
e for /i = 2/3. The behaviour of these curves is read- 
ily understood from the physics of Compton scattering. 
At low energies, the cross section is approximately that of 
Thomson scattering and a negligible fraction of the pho- 
ton's energy is transferred to the Compton electron. As a 
result, most of the incident flux is simply reflected back 
into space, with only a small fraction deposited within the 
reversing layer. 
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Fig. 1. Deposition of beam flux. The fractions of incident 
flux deposited above and below the photosphere (r = 2/3) 
are plotted against the photon energy e of the mono- 
energetic beam. The incident direction cosine /xo = 2/3. 

At high energies, on the other hand, the greatly re- 
duced Klein- Nishina cross section allows a large fraction 
of the incident photons to pass right through the reversing 
layer, so that most of the incident flux is thermalized be- 
low the photosphere. The local T e // is therefore increased 
and, in consequence, so also are temperatures in the re- 
versing layer. In this way, a sufficiently strong beam of 
7-rays with eo{MeV) <; 30 has an indirect effect on the 
observed line spectrum, with the n — 2 levels of He be- 
ing populated by the photosphcric continuum and by col- 
lisions with thermal electrons. This is the circumstance 
envisaged by Ferraro et al. (2003). 

However, for a broad range of intermediate photon en- 
ergies 0.04 <; eo(MeV) Ss 10, we see from the bell-shaped 
curve in Fig. 1 that a significant fraction of the incident 
flux is deposited in the reversing layer (r < 2/3) and thus 
directly affects line formation through impact excitations 
and ionizations. 

3.2. Curves of growth 

Fig. 1 shows that non-thermal excitation of the He I lines 
is favoured when the beam's photon energy is ~ 1 Mev. 
Accordingly, we now set eo = 1 MeV, since a definitive 
failure at this eo would be fatal for this interpretation. 

With e fixed, the single remaining parameter is 77, the 
conversion efficiency. We now take 77 to be a free parame- 
ter, postponing till later a discussion of observational and 
theoretical constraints on its magnitude (Sect. 4.2). 

As r\ increases from zero, the impact rates Tj in Eq.(5) 
increase proportionately - see Eq.(4), resulting eventually 
in column densities for the n = 2 levels sufficient for de- 
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Log 7) 



Fig. 2. Curve of Growth for Hci A5876A. The dimension- 
less equivalent width W\/\ is plotted against the conver- 
sion efficiency rj for the indicated beam parameters. The 
dashed curve corresponds to line formation by coherent 
scattering. 



tectable absorption of the photosphcric continuum. The 
simplest way of illustrating this effect is by constructing 
curves of growth with rj as the abscissa and dimcnsionlcss 
equivalent width W\/\ as ordinate. Accordingly, for Log q 
ranging from -5 to -0.5, the coupled radiative transfer and 
statistcal equilibrium problems of Sect. 2 have been solved 
and the resulting line profiles and equivalent widths W\ 
computed for all transitions represented in the adopted 
atomic model for He . 

Fig. 2 is the curve of growth for the line He I A5876A. 
We see that W\ increases linearly with rj for r\ ^ 4 x 10~ 4 
and thereafter more slowly. This is reminiscent of conven- 
tional curves of growth, where the linear part is followed 
by a flat part as the lines' Doppler cores become saturated 
while their damping wings remain optically thin. But here 
there is the additional effect of the partial filling-in of the 
absorption line by intrinsic line emission within the re- 
versing layer. These additional line photons are created 
in the radiative cascades following impact ionizations and 
excitations. 

The role of intrinsic emission in reducing W\ can 
be quantified approximately by also computing, for the 
given level populations, the line profile under the assump- 
tion of coherent scattering. For the adopted Schuster- 
Schwarzschild reversing layer, the residual intensity is then 
r„ = 1/(1 + 0.75Nj<iv), where Nj is the column density 
of absorbing atoms and a v is given by Eq.(6) . The curve 
of growth thus obtained is the dashed curve in Fig. 2. The 
vertical difference between the two curves increases with 
T], reflecting the growing contribution from cascades. 
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Fig. 3. Line profile of Hei A5876A for conversion efficiency 
77 = 0.01. The beam parameters are same as for Fig. 2. The 
dashed curve corresponds to line formation by coherent 
scattering. 

The actual line profiles used to derive these curves of 
growth are plotted in Fig. 3 for 77 = 10~ 2 . In this case, 
intrinsic emission reduces W\ from 0.35 to 0.28A. 

Since the strengths of other optical He 1 lines might 
also be relevant to testing this excitation mechanism, a 
selection of their curves of growth are plotted in Fig. 4. 
As with the A5876Aline, all lines have curves of growth 
that eventually depart from linearity in reponse to the 
combined effects of core saturation and intrinsic emission. 
The most extreme case is that of the A6678A line, whose 
W\ increases to a maximum of ~ 0.022A at 77 ~ 5 x 10~ 3 
and then declines precipitously, going into emission for 
T] > 1.3 x 10~ 2 . 

4. Comparison with observations 

In this section, the predictions of the non-thermal excita- 
tion model are compared with observational data for COM 
J1740-5340. 

4.1. The equivalent width of X5876A 

The prediction that each line has a maximum W\ offers an 
immediate prospect of observational contradiction. From 
Fig. 2, we see that reaches a maximum of 0.30A 

at rj = 4 x 10~ 2 . For comparison, observed values can be 
estimated from the spectra published by Sabbi et al. (2003; 
Fig. 4). These spectra show that W\ varies with phase but 
reaches strengths of ~ 0.19A at phases 0.02 and 0.56. 

This comparison shows that the predicted maximum 
exceeds observed values with a comfortable safety margin. 
This is just as well since, at these phases, a significant 
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Fig. 4. Curves of Growth for selected He I lines. The di- 
mensionless equivalent widths W\/\ are plotted against 
the conversion efficiency r\. The beam parameters are same 
as for Fig. 2. Note that the line A6678A is in emission for 
r) > 0.013. 

fraction of the companion's projected disk is not irradiated 
by the beam and is thus not expected to contribute to the 
integrated line profile. 

Eventually, when higher S/N spectra are obtained, it 
will be worthwhile to calculate line profiles integrated over 
the disk. This will not only give more definitive tests of 
line strengths but also allow the phase-dependent velocity 
displacement of the line centre from the orbital velocity of 
COM J1740-5340 to be tested. Such displacements are a 
consequence of only one hemisphere of the synchronously 
rotating companion (Ferraro et al. 2003) being irradiated. 
Of course, with data of exceptional quality, the technique 
of Dopplcr tomography could be used to make inferences 
about the pulsar's beam pattern. 

4.2. Conversion efficiency 

From Fig. 2, we find that IF5876 is comparable with the 
observed strengths for 77 <; 2 x 10 -3 . But this assumes the 
beam sweeps out solid angle = 1 steradian (Sect. 2.2). 
This must be increased by a factor <; 3.2 for the entire fac- 
ing hemisphere to be irradiated. Thus, more realistically, 
the required efficiency is rj <; 6 x 10~ 3 . The next question, 
therefore, is whether this is achievable for photon energies 
- 1 MeV. 

No data relevant to this question has yet been acquired 
for PSR J1740-5340. But such data is available for the 
binary millisecond pulsar PSR J0218+4232, whose spin 
period (2.3ms) and spin-down luminosity (2.5 x 10 35 erg 
s _1 ) (Mineo et al. 2000) are not too different from those 
of PSR J1740-5340. Detections and upper limits for pho- 



ton energies from 0.05 - 10,000 MeV are available for this 
pulsar. From this data, Kuiper et al. (2000) suggest that 
this pulsar's maximum luminosity density is reached in 
the COMPTEL MeV range just below the reported upper 
limits. Accepting this suggestion and thus interpreting the 
2cr upper limit reported for the COMPTEL 0.75-3 MeV 
band as a detection, we find rj b s ~ 15% for a 1 steradian 
beam. This is to be compared with rj b s ~ 7% found by 
Kuiper et al. from EGRET detections for energies between 
100 MeV and 1 GeV. 

This interpretation of the data for PSR J0218+4232 in- 
dicates that the efficiency required if eo ~1 MeV presents 
no problem. But note that, for this same data set, Dyks 
& Rudak (2002) offer a theoretical interpretation which 
predicts fluxes in the COMPTEL bands far below the re- 
ported upper limits. Specifically, at ~1 MeV their model 
predicts r\ ~ 10~ 3 , corresponding to a barely detectable 
line at A5876A. 

4.3. Heating 

The 1-steradian efficiency n ~ 2 x 10~ 3 needed to produce 
the observed strength of the A5876A line corresponds to 
an incident 7-ray flux Tq = 7.5 x 10 8 erg cm~ 2 s _1 . Even 
if all of this were thcrmalizcd below the photosphere, the 
increase in T e // of the irradiated hemisphere is only ~ 
20K. Thus the proposed explanation of the A5876A line 
is consistent with the absence of a detectable photometric 
signature of heating (Sect. 1). 

4.4. Hydrogen lines 

Hydrogen line profiles have also been calculated as a 
function of 77. The initial thought was that the observed 
strength of the A5876A line might be achieved simulta- 
neously with Ha going into emission, thus explaining the 
star's earlier listing as a BY Dra candidate (Taylor et al. 
2001). In the event, at values of 77 corresponding to A5876A 
reaching its observed strength, the contribution of intrin- 
sic emission to Ha is still neglible. 

However, this is not a failure of the model since it is 
clear from the spectra published by Sabbi et al. (2003) 
that Ha emission originating close to the surface of COM 
J1740-5340 is not due to irradiation. This conclusion fol- 
lows because the Ha emission remains stong even when 
most of the facing hemisphere is hidden from the observer. 
Most likely, as Sabbi et al. conclude, the emission has a 
chromospheric origin as commonly found in short-period, 
low mass binaries. 

4.5. The equivalent width of X6678A 

The above discussion has focussed on He 1 A5876A since 
this the line from which Ferraro et al. (2003) and Sabbi 
et al. (2003) infer that the pulsar's beam is highly 
anisotropic. 
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Of the other He I lines, the one whose predicted 
strength is most likely to be in conflict with observation 
is the A6678A line, which is the analogue for the sin- 
glets of the triplet line A5876A. Ferraro et al. (2003) re- 
port the presence of the A6678A line but not its strength. 
According to these calculations, its maximum equivalent 
width is only 0.022A, which is probably inconsistent with 
detetection. 

Subsequent to the original submission of this pa- 
per, E.Sabbi (Bologna) has confirmed that the A6678A 
line is indeed reliably detected. Moreover, it reaches an 
equivalent width well in excess of the above maximum. 
Specifically, at phase 0.56, she reports that 1^6678 — 0.19 
A, almost an order of magnitude larger than the predicted 
maximum. 

The weakness of A6678A in these calculations contrasts 
with its strength in models for Type lb SNe (Paper I, 
Fig.l). The origin of this difference is chemical composi- 
tion. The atmospheres of Type lb SNe are hydrogen-free, 
so there is no loss of population from singlet levels due to 
the ionization of hydrogen by photons emitted by the per- 
mitted ground state transitions n 1 P° — > l x 5 (Sect. 2.7). 

This explanation for the predicted weakness of the 
A6678A line is confirmed by recomputing the model at 
Logrj = —2.3 but now with the line escape probabilities 
of Sect. 2.7 = 0. The result is a dramatic increase in 
W-Wg from 0.022 to 0.33A. 

Interestingly, in this experiment W5876 also increases, 
from 0.26 to 0.35A. The reason for this is the removal 
of the population decrement of the n = 2 singlet levels 
relative to the corresponding triplet levels. In consequence, 
there is no longer a significant net collisional transfer of 
excitation from the triplets to the singlets. 

4. 6. Other lines 

In his report, the referee F.R. Ferraro (Bologna) generously 
included unpublished information on other He lines. This 
data, as also that from E. Sabbi discussed above, is from 
the spectra described by Ferraro et al. (2003) and Sabbi 
et al. (2003). 

The triplet line A3889A has a predicted strength close 
to that of the A5876A line - see Fig. 4 - and so should be 
seen were it not blended with the H8 Balmer line. 

The triplet line A4471A has a predicted equivalent 
width of 0.06A when W^s76 = 0.19A, as observed at phases 
0.02 and 0.56. Though weak, a line of this strength should 
still be visible, but Ferraro and colleagues finds no evi- 
dence of it. This is undoubtedly a difficulty but one un- 
likely to be specific to the non-thermal mechanism. The 
reason for this is that the A5876 and A447lA lines have the 
same lower level and the ratio of their predicted strengths 
is approximately equal to the ratio of their f- values. Thus 
a similar ratio should be expected for the thermal mecha- 
nism proposed by Ferraro et al. (2003) 



5. Conclusion 

The aim of this paper has been to investigate the possi- 
bility that the formation mechanism for He 1 absorption 
lines in the spectrum of COM J1740-5340 is similar to that 
operating in the atmospheres of Type lb SNe. This alter- 
native merits serious consideration in view of the powerful 
and remarkable constraint on the pulsar's radiation pat- 
tern that follows if the He I lines are indeed formed by 
conventional thermal processes. 

The result of this investigation is moderately discour- 
aging for the non-thermal interpretation but not yet a 
decisive rejection. On the positive side, the crucial line 
A5876A achieves the observed strength at seemingly fea- 
sible values of the conversion efficiency r\. Moreover, the 
corresponding heating effect is consistent with the non- 
detection of any such effect photometrically. 

However, despite these successes, the strengths of other 
Hei lines, especially A6678A, present difficulties. The ques- 
tion then is whether such difficulties are fatal or can 
be overcome by refinement or modification of the simple 
model used here. Among refinements to be considered are 
the effect of the beam on the continuum opacity, specifi- 
cally the dissociation of the H~ ion by impacts (Geltman 
1960), the ionization of surface layers by impacts, and the 
inclusion of an empirical model for the chromosphere dis- 
covered by Sabbi et al. (2003). A higher level treatment of 
a stellar atmosphere's response to irradiation by gamma 
rays is an attractive technical challenge and likely to have 
application to other close companions of exotic objects. 

Nevertheless, given the relatively low S/N of the spec- 
tra of COM J1740-5340 and the poorly constrained pa- 
rameters of the irradiation model, it is perhaps doubt- 
ful that a decision between the thermal and non-thermal 
line formation mechanisms can soon be reached purely 
on the basis of line diagnostics. Accordingly, the predic- 
tion by Sabbi et al. (2003) that up to 1% of the disk of 
COM J1740-5340 radiates with T eff ~ 10,000K should 
be tested with deep UV imaging at A 2000A. A de- 
tection of a photometric signature of a heating effect at 
such wavelengths would be a decisive confirmation of the 
thermal mechanism and therefore of a narrowly-confined 
pulsar beam. 

Acknowledgements. I am indebted to E.Sabbi and the referee, 
F.R. Ferraro, for unpublished data on He line strengths. 
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